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Abstract 


Hydrogen production via steam electrolysis may involve less electrical energy consumption than conventional low temperature water electrolysis, 
reflecting the improved thermodynamics and kinetics at elevated temperatures. The present paper reports on the development of a one-dimensional 
dynamic model of a cathode-supported planar intermediate temperature solid oxide electrolysis cell (SOEC) stack. The model, which consists of 
an electrochemical model, a mass balance, and four energy balances, is here employed to study the steady state behaviour of an SOEC stack at 
different current densities and temperatures. The simulations found that activation overpotentials provide the largest contributions to irreversible 
losses while concentration overpotentials remained negligible throughout the stack. For an average current density of 7000 A m~? and an inlet 
steam temperature of 1023 K, the predicted electrical energy consumption of the stack is around 3 kW h per normal m° of hydrogen, significantly 
smaller than those of low temperature stacks commercially available today. However, the dependence of the stack temperature distribution on the 
average current density calls for strict temperature control, especially during dynamic operation. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Hydrogen is regarded as a leading candidate for alternative 
future fuels. It has the potential to address the environmental and 
energy security issues associated with fossil-derived hydrocar- 
bon fuels. Among many hydrogen production methods, water 
electrolysis is a well-established technique, which is capable of 
producing emissions-free hydrogen if used in conjunction with 
renewable or nuclear energy. However, water electrolysis has 
not had a significant commercial penetration, mainly due to its 
high electricity consumption and thus high operating or running 
cost [1]. 

Steam electrolysis at elevated temperatures might offer a 
solution by consuming less electricity than is required at ambi- 
ent conditions. The minimum energy consumption during steam 
electrolysis is evaluated as a function of temperature in Fig. 1, 
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using thermodynamic tables [2]. Here the gas compositions are 
taken to be 50 mol% H2/50mol% H20 and 100 mol% Odo. In 
the figure, TAS and AG represent, respectively, the thermal and 
electrical energies consumed in the reaction, while AH cor- 
responds to the total energy consumption, in the absence of 
irreversible losses. Although AH remains approximately con- 
stant, AG tends to decrease and TAS tends to increase as the 
temperature is raised, providing more opportunities to apply 
waste heat from other processes, or to make effective use of the 
heat generated via irreversible losses. At 1273 K, almost 30% of 
the total energy consumption can be satisfied through the sup- 
ply of thermal energy alone. The increase in temperature can 
also eliminate the need for expensive catalysts, which may be 
required for some low temperature water electrolysers [3]. How- 
ever, in the design of a solid oxide electrolysis cell (SOEC), there 
are other factors that need to be taken into account, which influ- 
ence the performance of the actual device as will be explained 
in the following sections. 

Steam electrolysis is performed using an SOEC, which can be 
seen in simple terms as the reverse operation of a solid oxide fuel 
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Nomencaure Te temperature of the cathode gas stream at the inlet 
Cy,, Cmo concentration of hydrogen and steam in the (K) 
cathode gas stream (mol m7?) uc, ua velocity of the cathode and anode gas streams 
Ci concentration of gas species i in the cathode and (ms!) 
anode gas streams (mol m~?) U cell potential (V) 


Cp,C», Cpa» Cp,s, Cpr heat capacity of the cathode and 
anode gas streams, solid structure and intercon- 
nect (J kg7! K7!) 


Cc? concentration of gas species i at the cell inlet 
(mol m~?) 
Cees CHO concentration of hydrogen and steam at the 


triple phase boundary (mol m~?) 

dnc, dna hydraulic diameter of the cathode and anode gas 
channels (m) 

Deffcathode average effective diffusivity coefficient of the 
cathode (m? s7!) 

Ecathode, Eanode activation energy for the exchange current 
density of the cathode and anode (J mol`!) 

Electrode activation energy for the exchange current den- 
sity of an electrode (J mol`!) 


F Faraday’s constant (C mol`!) 
AG Gibbs free energy change of the reaction (J mol~! 
H2) 


hc, ha, hs, hy height of the cathode and anode channels, 
and thickness of the solid structure and intercon- 


nect (m) 
AH enthalpy change of the reaction (J mol~! H2) 
J local current density (A m7?) 


Jo,cathode> jO,anode eXchange current density of the cathode 
and anode (A m~?) 

Jo,electrode exchange current density of an electrode 
(A m~?) 

kc, ką convective heat transfer coefficient between the 
solid parts of the cell and the cathode and anode 
gas streams (J m~? s7! K7!) 

kcathode, Kanode pre-exponential factor for the exchange 
current density of the cathode and anode 
(Q7! m~?) 

Kelectrode pre-exponential factor for the exchange current 
density of an electrode (Q7! m?) 


L cell length (m) 

Nuc, Nua Nusselt number of the cathode and anode gas 
streams 

Po, oxygen pressure in the anode gas stream (Pa) 

PYP oxygen pressure at the triple phase boundary (Pa) 

R reaction rate (mol si m~?) 


Ronm _ total resistance of the cell including electric and 
ionic resistances (Q m?) 

R gas constant (J mol`! K7!) 

AS entropy change of the reaction (J K~! mol! H2) 

t time (s) 

T temperature (K) 


y? standard potential (V) 

yy reversible potential (V) 

urev.TPB reversible potential estimated for the gas concen- 
trations at the triple phase boundary (V) 

W cell width (m) 

x axial coordinate (m) 


Greek letters 

a transfer coefficient 

€s,&€] emissivity of the solid structure and interconnect 

Nact,cathode> Nact,anode activation overpotential losses at the 
cathode and anode (V) 

Nconc total concentration overpotential losses (V) 

Neonc,cathode> Nconc,anode Concentration overpotential losses 
at the cathode and anode (V) 

NOhm Ohmic losses (V) 

Ntotal total irreversible losses (V) 

Ac, ÀA, As, Ay thermal conductivity of the cathode and 
anode gas streams, solid structure and intercon- 
nect (J m~! s7! K7!) 

vi stoichiometric coefficient of gas species i 

PC, PA, Ps, pr density of the cathode and anode gas 
streams, solid structure and interconnect (kg m~?) 

o Stefan-Boltzmann constant (W m~? K~*) 

Ocathode; Sanode Electric conductivity of the cathode and 
anode (Q7! m7!) 

Gelectrolyte ionic conductivity of the electrolyte (Q27! m7 1) 

Tcathodes Telectrolytes Tanode thickness of the cathode, elec- 
trolyte and anode (m) 


Energy consumption (kJ.mol! H2) 


373 573 773 973 1173 1373 
Temperature (K) 


Fig. 1. Thermal (TAS), electrical (AG) and total (AH) energy consumption 
during steam electrolysis as a function of temperature. The gas compositions 
are taken to be 50 mol% H2/50 mol% H20 and 100 mol% Oz. Irreversible losses 
are not taken into account. 


Tc, Ta, Ts, Tj temperature of the cathode and anode gas 
streams, solid structure and interconnect (K) 
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Table 1 
Chemical reactions involved in an SOEC 


Cathode H20 +2e7 —> H2 +027 
Anode o% > 30» +207 
Overall H.,O > H2 + 302 


cell, allowing for the opportunity to apply recent developments 
in SOFCs to the field. An SOEC consists of a three-layer solid 
structure (composed of porous cathode, electrolyte and porous 
anode) and an interconnect plate. Steam is introduced at the cath- 
ode side of the solid structure where it is reduced into hydrogen, 
releasing oxide ions in the process. The oxide ions then migrate 
through the electrolyte to the anode where they combine to form 
oxygen molecules, releasing electrons. The balanced chemical 
equations for the reactions are shown in Table 1. Although the 
use of proton conductors has also been investigated [4—6], yttria- 
stabilised zirconia (YSZ), which is an oxide ion conductor, is 
generally used for the electrolyte in SOECs. Typical materials 
for the cathode are nickel-YSZ cermets and those for the anode 
are perovskite oxides such as lanthanum manganite. 

In SOFC technologies, there is increased interest in interme- 
diate temperature SOFCs (IT-SOFCs), which typically operate 
between 823 and 1073 K, to allow for a wider range of mate- 
rials, more cost effective SOFC fabrication [7] and increased 
durability. The reduction in operating temperature has also been 
seen in the field of SOECs in the limited number of experimen- 
tal projects reported over the recent decades [8—10], where the 
temperature at which the cells are tested has decreased from 
1273 to 1073 K. Initially SOECs were developed as tubular cell 
structures to avoid sealing problems, which is a major issue in 
segregating the produced H2 and O2 in planar cells, particularly 
over multiple thermal cycles [11]. Today, despite such sealing 
difficulties, the trend in cell design is to employ planar structures. 
A planar structure permits high packing density and significantly 
smaller hot volume in the system than that allowed by a tubular 
design [10]. The main motivation behind the use of planar cells 
is, however, the low manufacturing cost [12]. The planar design 
also allows for shorter current path, reducing the Ohmic resis- 
tance within the cell [12]. A planar SOEC can be categorised 
as being either an electrolyte-supported or electrode-supported 
design. The former employs the electrolyte as the support struc- 
ture and is suitable for high temperature operation in which the, 
often large, Ohmic resistance associated with a thick electrolyte 
can be reduced. In an electrode-supported cell, on the other hand, 
one of the electrodes is the thickest part of the solid structure. 
Such a design has been developed in order to minimise Ohmic 
resistances in SOFCs operating at intermediate temperatures 
and may also be applied to intermediate temperature SOECs 
(IT-SOECs). 

A mathematical model is an important design tool for such 
devices as an SOEC, which are still in the development stage. It 
allows the prediction of the behaviour of the device under dif- 
ferent process conditions and assists in the optimisation of its 
performance. In particular, such a model is essential in under- 
standing the response of an electrolyser under steady electrical 
power input, such as those from the grid or nuclear energy, as 


well as under an intermittent renewable electrical power input. 
The present paper reports the development of a dynamic model 
of a conceptual IT-SOEC produced within the gPROMS mod- 
elling environment [13]. The paper introduces the IT-SOEC 
model and the results describing the simulated steady state 
behaviour of such an electrolyser. Issues related to thermal man- 
agement of the cell are discussed. In spite of the model’s dynamic 
capability, only the steady state results are discussed here. The 
dynamic response will be presented in a subsequent publication. 


2. IT-SOEC mathematical model 


Today, although studies on the thermodynamics of electroly- 
sis at elevated temperatures can be found in several publications 
[9,14,15], mathematical modelling of an SOEC stack is still not 
an active area of research. To the authors’ knowledge, neither 
results related to the spatially distributed stack performance nor 
the transient behaviour of a stack have been published. Here, 
a previously developed dynamic model of an anode-supported 
planar IT-SOFC stack [7,16] has been modified to produce a 
one-dimensional model for a planar IT-SOEC stack that oper- 
ates across a wide range of temperatures, using a thick film YSZ 
electrolyte in a cathode-supported cell configuration. 

To ensure a sufficient rate of Hy production, an SOEC sys- 
tem must consist of several repeating cells assembled in stacks. 
However, the models of such stacks are usually constructed for 
the smallest unit cell, which is assumed to describe the response 
of the whole stack subject to the use of adequate boundary con- 
ditions. Here, the modelled unit cell is considered to be in the 
centre of a large stack such that there are no end effects. Although 
interconnects normally provide the gas flow channels, above and 
below the solid structure, the effect of individual passages is 
here neglected. The pressure drop along the gas channels is also 
assumed negligible at the operating pressure of 0.1 MPa. For the 
purpose of modelling, the unit cell is considered to be composed 
of four components, the cathode and anode gas streams, solid 
structure and interconnect. Fig. 2 shows the schematic view of 
such a unit cell. The model developed consists of an electro- 
chemical model, a mass balance for the cathode stream, and 
four energy balances for the cathode and anode streams, solid 


H, / H,O —> 
[ae 


Fig. 2. Schematic view of a planar SOEC stack. 


— H,/H,0 
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structure and interconnect. In spite of the temperature gradi- 
ent expected along the cell, the properties of the gas streams, 
solid structure and interconnect are assumed to be constant in 
order to achieve a more computationally tractable model. This 
assumption is reported valid for the previously developed IT- 
SOFC model during steady state simulations [17] and is thus also 
considered applicable for the IT-SOEC model. Finally, ideal gas 
behaviour is assumed for the cathode and anode streams where 
appropriate. 


2.1. Electrochemical model 


An electrochemical model is used to relate such variables 
as gas species concentrations, cell component temperatures and 
average current density to the electrical potential of the cell, 
which can then be used to calculate the electrical energy con- 
sumption. The electrical potential is taken as constant along the 
cell as the electrodes are assumed to be good electric conduc- 
tors, and thus it is not a local quantity. As shown in Eq. (1), the 
cell potential corresponds to the sum of the reversible potential 
and all the irreversible losses that occur as the electrical cur- 
rent is passed through the cell. Such irreversible losses include 
Ohmic losses, concentration overpotentials and activation over- 
potentials, which are all partly responsible for the heat produced 
within the cell. In general, irreversible losses are dominated at 
high current densities by concentration overpotentials and at low 
current densities by activation overpotentials. 


U= U(x) + Nohm(x) + Neonc,cathode(X) + Neonc,anode(*) 
+ Nact,cathode (x) + Nact,anode (x) (1) 


The reversible potential is determined through the Nernst 
equation, which predicts the minimum electrical potential 
required to split H20 at a particular location along the cell with 
a specific temperature and gas concentrations. Eq. (2) has been 
derived assuming ideal gas behaviour. The solid structure tem- 
perature is applied in the Nernst equation as the reactions take 
place at the triple phase boundaries (TPBs) where the electrolyte, 
porous electrode and gas species meet at the interfaces between 
the electrolyte and the electrodes. U? represents the standard 
potential, which corresponds to the reversible potential calcu- 
lated at unity activity. It is determined as a function of the Gibbs 
free energy change involved in the reaction, taking the partial 
pressures of H2, H20 and O; to be at standard pressure. 


—5,0.5 
RTs) | | Cra CLPo, x 10 5] 


rev — 770 
U(x) = U (x) + ZF E 


(2) 


2.1.1. Ohmic losses 

Ohmic losses originate from the resistance to the conduc- 
tion of ions through the electrolyte and electrons through the 
electrodes and interconnect, as well as from the contact resis- 
tances, which exist between cell components. Ohmic losses are 
linearly proportional to the current density and can be deter- 
mined through Ohm’s law as shown in Eq. (3), where j and 
Ronm represent the local current density and local cell resistance, 


respectively. Ronm is calculated using the conductivity data of 
each layer of the solid structure. In Eq. (4), cross plane charge 
flow and series connection of resistances have been assumed 
[7]. The resistance through the interconnect and the contact 
resistances are considered negligible. The conductivity of the 
electrodes are assumed constant along the cell while that of the 
electrolyte is estimated as a function of local temperature [18]. 


Nom) = JŒ) Ronm) (3) 
Tcathode Telectrolyte Tanode 

Rohm(*) = (4) 
Ocathode Celectrolyte (x) Oanode 


2.1.2. Concentration overpotential losses 

Concentration overpotentials occur if mass transport of the 
gas species involved in the reaction becomes the limiting factor. 
When current is flowing in the cell, reactant H20 is consumed 
and product H3 and O3 are formed at the TPBs. If the fluxes of the 
reactant to, and the products away from, the TPBs are slower than 
that corresponding to the discharged electrical current, concen- 
tration gradients for the species develop in the porous electrode. 
The concentration overpotentials are the result of such concen- 
tration gradients causing the actual reversible potential at the 
TPBs, shown in Eq. (5), to deviate from the reversible poten- 
tial predicted by the Nernst Eq. (2) using bulk concentrations in 
the cathode and anode streams. High current densities and high 
steam utilisation factors can contribute in creating significant 
concentration overpotentials. 


—5,0.5 
YTB (x) — Gy 4 RPS, | Cr LPO,” x 1071 
? x)= X 
2F CHB (x) 
(5) 
w = BRO, Ci (CHOC) 
er ae Cry (x) CHP (x) 
RTs), | POPP 
In |2 6 
+ GF" | Po, = 


The difference between Eqs. (2) and (5) corresponds to the 
concentration overpotential (6) where the first and second terms 
on the right hand side represent the concentration overpotential 
at the cathode and anode, respectively. The cathode concen- 
tration overpotential is presented separately in Eq. (7) while the 
anode concentration overpotential is not considered here, assum- 
ing the difference between the O2 pressure at the TPBs and that 
in the bulk stream to be negligible. Eqs. (8) and (9) show the 
TPB concentrations, determined for equimolar counter-current 
one-dimensional diffusion of H2 and H20. Dett.cathode represents 
the average effective diffusivity coefficient in the cathode, con- 
sidering a binary gas mixture of Hz and H20. A method for 
expressing the mass diffusion activities in an SOFC [19,20] has 
been adopted for an SOEC in deriving Eqs. (8) and (9). 


RTO CE’ œC oa) 7 
Nconc, cathode (x) = ORF n Ca) CIPB (x) (7) 
2 2 
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Table 2 
Pre-exponential factor and activation energy for the exchange current density 


Keathode 654 x 10° Q7! m~? E cathode 140 x 103 Jmol7! 
Kanode 235 x 10° Q7! m~? Emnode 137 x 103 Jmol”! 
TPB Tcathode . 
Cy, œ) = Cu, (x) + > j(x) (8) 
= 2F 'Deff,cathode 
CIPB Tcathode ; 
Cho) = Cmo) — z ~~~ _ jx) (9) 
2 2F. 'Deff,cathode 


2.1.3. Activation overpotential losses 

Activation overpotentials are related to the chemical kinet- 
ics of the reactions. They are classically determined through 
the Butler-Volmer equation. For the cathode activation over- 
potential, the extended form of the Butler-Volmer equation is 
applied in Eq. (10) to account for the differences in the gas con- 
centrations at the TPBs and in the bulk stream. For the anode 
activation overpotential (11), the simplified form is used, as the 
difference between the O2 concentration at TPBs and that in the 
bulk stream has been assumed negligible. Although reactions are 
generally rapid and activation overpotentials tend to be small at 
high temperatures, they can become the most significant form 
of irreversible losses as the operating temperature is reduced. 


es - CHEB (x) 2(1 — a) F a 
JX) = J0,cathodelX Ca x ex | NTs) Nact,cathode »| 
CIPB 
Wa —2aF 
Cho 
Croa) JP Ts rata) | (10) 
PE 21-—a)F 
Jx) = jo,anode(x) | exp Ra enod) 
—2aF 11 
p E rasan] (11) 


a (usually taken to be 0.5), jocathode and jo,anode represent 
the transfer coefficient and the exchange current density for 
the cathode and anode, respectively. Electrode exchange current 
densities are expressed in Eq. (12) as a function of the pre- 
exponential factor and activation energy. Here it was assumed 
that in an SOEC, the forward and backward reactions at each 
electrode as well as the electrode materials are the reverse of 
those in an SOFC. Hence the values used for Kanode and Eanode 
in previous IT-SOFC studies have been substituted into Keathode 
and Ecathode in the IT-SOEC model. The same assumption is used 
when assigning the values to the IT-SOEC anode kinetic param- 
eters. The parameters employed for the IT-SOEC simulations 
are shown in Table 2. 


HT s(x) 
2F 
electrode € {cathode, anode} 


J0,electrode (x) = 


RTs) 


—E electrode 
(12) 


Kelectrode €Xp | 


2.2. Mass balance 


The composition of the cathode stream evolves along the 
cell as the H2 /H20 mixture moves towards the outlet. Such a 
composition change is tracked in the model using the cathode 
stream mass balance (13) and the boundary conditions at the cell 
inlet (14), predicting the H? and H20 concentrations at each 
location along the channel. R represents the rate of reaction, 
which is related to the local current density through Faraday’s 
law (15), assuming a current efficiency of 100%. As for the anode 
channel, it is assumed that any O2 produced flows naturally out 
of the channel and thus there is no forced gas movement as in the 
cathode side. Therefore, the developed model does not consider 
an anode mass balance. Instead, the anode stream is modelled 
as a flow of pure O2 at constant pressure along the cell. The 
cathode and anode stream velocities are assumed constant and 
determined from the outlet flow rate of the streams. 


a a 1 
LO! = ucla) + pma i € {H2, H20} 
(13) 
Ci(0) = C?, ie{H2, H20} (14) 
R(x) = Mm (15) 


2.3. Energy balances 


Eqs. (16)—(23) are used in predicting the temperatures of the 
cathode and anode streams, solid structure and interconnect. In 
Eqs. (16) and (18), the first term on the right hand side relates to 
the energy carried in the gas streams as they flow along the chan- 
nels. The energy balances consider the thermal fluxes between 
the gas streams and the solid parts of the cell to be fully described 
by convection. The convective heat transfer coefficients for the 
cathode and anode streams (kc and ka) are calculated using Eq. 
(24), assuming constant Nusselt number and hydraulic diameter 
(25). In the cathode stream, kc is assumed equivalent for both 
the convective heat transfer over the solid structure and intercon- 
nect surfaces. The same assumption is applied to ka in the anode 
stream. The thermal fluxes along the solid parts of the cell are 
modelled using Fourier’s law of heat conduction while radiative 
heat exchange is taken into account between the solid structure 
and interconnect. The entire enthalpy change of the reaction is 
assumed to occur in the solid structure. Mean density and mean 
heat capacity values are applied in the anode and cathode energy 
balances. These are calculated as spatial averages of all the val- 
ues along the cell. The remaining gas and material properties 
are taken to be constant along the cell. 


ð 0 Ke 
gl Teo] = —uc—[Te(x)] + [Ts(x) — Te(x)] 
t ox Pcep,chc 
kc 
+ [Tr(x) — Te(x)] (16) 
Pcep,chc 


Tc(0) = TE (17) 


132 J. Udagawa et al. / Journal of Power Sources 166 (2007) 127-136 


a a ka 
zl TA] = —ua—[Ta(x)] + [Ts(x) — Ta] 
t Ox PACp, ANA 
ka 
[T — Ta(x)] (18) 
PACp ANA 
if T 
Ta(0) = BOTAO (19) 
Ts(x)) = 2S a T. keg, T, 
zl s(x)] = ae zl s(x)] rar s(x) — Tce(x)] 
ka 
[Ts(x) — Ta(x)] 
PsCp,shs 
2 o[Ts(x)* — T(x)*] 
Pscp,shs | (1/es) + (1/21) — 1 
———— [—AA (x) R(x) + j(x)U] (20) 
PSCp,shs 
a a 
9, sl = 0, —[Ts(L)] = 0 (21) 
x ox 
meas : imo- “81 — Te] 
a x = x Xx X 
a | pICp,1 Ox? i PrCp thy i S 
[Ty(x) — Ta(x)] 
PIC Ar 
2 o[Ts(x)* — T(x)*] (22) 
prcp,rhy | (1/es) + (1/21) — 1 
a a 
—[T%(0)] = 0, —[N(L)] = 0 (23) 
ox Ox 
Wye — Nag A 
kc = Nuc he ka = Nuh (24) 
= 2Whc _ 2Wha 
dnc = Wine’ bA = Woh, F (25) 


3. Simulation results and discussion 


The system of partial differential and algebraic Eqs. (1)-(4) 
and (7)—(25) presented in the previous sections is solved via the 
finite difference method using gPROMS Model Builder 2.3.5 
[13]. The model input parameters and the operating conditions 
applied for the steady state simulations are introduced in Sec- 
tion 3.1. Section 3.2 discusses the effects of current density and 
temperature on the cell potential while Section 3.3 presents the 
spatially distributed simulation results. 


3.1. Input parameters and inlet operating conditions 
In addition to the input parameters such as cell geometry, 


material properties and inlet stream conditions, the IT-SOEC 
model described in the previous sections involves two degrees 


Table 3 

Model input parameters and operating conditions 

Cathode channel height, hc 0.001 m 
Anode channel height, ha 0.001 m 

Solid structure thickness, hs 570 x 107° m 
Interconnect thickness, hy 500 x 107° m 
Cell length, L 0.4m 

Cell width, W 0.1m 
Cathode thickness, Tcathode 500 x 1076 m 
Electrolyte thickness, Telectrolyte 20 x 107° m 
Anode thickness, Tanode 50 x 107 m 
Cathode electric conductivity, Ocathode 80 x 10° Q7! m7! 


33.4 x 10° exp 

(—10.3x 10°/Ts) Q7! m7! 
8.4 x 10° Q7! m7! 

36.6 x 1076 m~? s7! 


Electrolyte ionic conductivity, Celectrolyte 


Anode electric conductivity, O anode 
Cathode average effective diffusivity, Deff,cathode 


Solid structure emissivity, £s 0.8 

Interconnect emissivity, £1 0.1 

Solid structure heat capacity, Cp,s 500Jkg-! K7! 
Interconnect heat capacity, Cp 1 500Jkg-! K7! 
Solid structure thermal conductivity, As 2Jm7!s7! K7! 
Interconnect thermal conductivity, Ay 25Jm7! s7! K7! 
Solid structure density, ps 5900 kg m~’ 
Interconnect density, py 8000 kg m~? 
Cathode stream Nusselt number, Nuc 3.09 

Anode stream Nusselt number, Nua 3.09 

Transfer coefficient, œ 0.5 

Cathode stream inlet temperature, Te 1023K 

Cathode stream inlet composition 10 mol% H2/90 mol% H20 
Anode stream composition 100 mol% O2 
Operating pressure 0.1 x 10° Pa 
Average current density 5000 A m~? 
Steam utilisation factor 80% 


of freedom. Here, average current density and steam utilisation 
factor are adopted as the user-defined parameters from which 
the inlet molar flux of steam is determined. The average cur- 
rent density is directly proportional to the H? production rate 
while the steam utilisation factor is the fraction of total inlet 
steam consumed by the reaction. Although a low steam utilisa- 
tion factor is preferred from the stack efficiency point of view, 
a reduction in the steam utilisation factor results in an increased 
cathode stream flow rate for a given H2 production rate. In the 
system design, the selection of the steam utilisation factor must 
take into account the efficiency of the entire system as well as the 
added investment cost involved in accommodating the increased 
cathode stream flow rate. It is also important for the steam utili- 
sation factor to be low enough to avoid a significant increase in 
the cathode concentration overpotential caused by steam starva- 
tion near the cell outlet. Here, a value of 80% is chosen. Table 3 
presents the model input parameters selected for the steady state 
simulations including the cell geometry, the material proper- 
ties, the cathode stream inlet temperature and composition, the 
operating pressure, the average current density and the steam 
utilisation factor. Due to the oxidising environment that pure 
steam would create at elevated temperatures, such a gas compo- 
sition is not recommended for the cathode stream in an SOEC. 
Instead, conceptual system designs may consider the recircu- 
lation of some of the product hydrogen, which is mixed with 
the steam at the inlet to ensure reducing conditions [10,14,21]. 
Here, 10 mol% Hp in the cathode inlet stream is assumed suf- 
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Fig. 3. Reversible potential and irreversible losses as a function of current 
density at 1023 K. The gas compositions are equivalent to those at the cell inlet. 


ficient [10]. Note that the model assumes the cell with a total 
length of 0.4 m, comprising of four individual parts connected in 
series, each with a dimension of 0.1 m x 0.1 m. The cell geom- 
etry and material properties are based on the values employed 
in a previously developed IT-SOFC model [7]. 


3.2. Effects of operating conditions on the cell potential 


The selection of operating conditions for an SOEC can 
significantly influence the irreversible losses, altering the cell 
potential required for the electrolysis. Such effect is illustrated 
in Figs. 3 and 4, respectively, showing the cell potential as a func- 
tion of current density and temperature. Note that these curves 
do not result from the full IT-SOEC model described in Section 
2. Instead, Figs. 3 and 4 are produced using only the electro- 
chemical model in Section 2.1, taking the cathode and anode 
stream compositions to be equivalent to those at the cell inlet. 

Fig. 3 shows the predicted cell potential at 1023 K, indicating 
the contributions of the reversible potential and each form of irre- 
versible losses. As can be seen, all the irreversible losses increase 
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Fig. 4. Reversible potential and irreversible losses as a function of temperature 
at 5000 A m~?. The gas compositions are equivalent to those at the cell inlet. 


with current density. At 8000 A m~?, the cell potential is around 
1.28 V with the anode activation overpotential being the largest 
form of irreversible losses at around 0.18 V, followed by the 
Ohmic losses and the cathode activation overpotential at around 
0.11 V and 0.09 V, respectively. In spite of the cathode-supported 
nature of the cell, the cathode concentration overpotential con- 
tributes less than 0.02 V to the cell potential at 8000 A m~?. 
While the Ohmic losses show essentially linear increase with the 
current density, the activation overpotentials follow non-linear 
paths. Although a decreased current density allows the reduc- 
tion in the cell potential and consequent decrease in electrical 
energy consumption, a low current density also results in a low 
Hy production rate per unit cell area. Therefore, the selection 
of the current density at which the cell should be operated must 
consider the operating cost involved with the electrical energy 
consumption of the stack as well as the investment cost associ- 
ated with the stack area required in supporting the specified H2 
production rate. 

Another operating condition to be selected in the design of 
an SOEC is the operating temperature. The reversible potential 
and irreversible losses are illustrated as a function of tempera- 
ture in Fig. 4. Except for the cathode concentration overpotential, 
which remains approximately constant, all the irreversible losses 
decrease with temperature. Between 973 and 1073 K, the cell 
potential is significantly decreased from around 1.36 to 1.03 V, 
with the combined irreversible losses being reduced by more 
than 60%. As expected from thermodynamics, the reversible 
potential given by the Nernst equation also decreases with tem- 
perature. The selection of the operating temperature must again 
consider the operating cost linked with the electrical energy 
consumption of the stack. In addition, it must account for 
the investment cost involved in creating both a stack with the 
required thermal stability as well as the whole system that needs 
to operate at that specified temperature. If the use of waste heat 
from external sources is considered, the temperature of such a 
heat stream would also become an important factor. 


3.3. IT-SOEC operation 


This section presents the spatially distributed steady state 
simulation results obtained by solving the full IT-SOEC model 
described in Section 2. Profiles of the temperature, the reversible 
potential, the irreversible losses and the local current density are 
illustrated along the cell. Different modes of operation for an 
SOEC are explained and the issues associated with the thermal 
management of the cells are discussed. 

Fig. 5 shows the temperature distribution along the cath- 
ode stream for the cell operation at average current densities of 
5000 A m~? and 7000 A m7, maintaining all other conditions 
consistent with Table 3. Note that the cathode stream tempera- 
ture changes rapidly after the cell entry, reaching the temperature 
of other cell components within a very short distance of the inlet. 
As can be seen in Fig. 5, the temperature gradient depends sig- 
nificantly on the average current density. While the temperature 
decreases along the cell at 5000 A m2, the heat accumulation 
causes the cell temperature to increase towards the outlet in the 
7000 A m~? case. This dependence of the cell temperature gra- 
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Fig. 5. Cathode stream temperature along the cell for average current densities 
of 5000 A m~? and 7000 A m~? and the conditions in Table 3. 


dient on the average current density, or the hydrogen production 
rate, thus calls for a strict cell temperature control. 

The effect of current density on the temperature gradient can 
be explained by considering the thermal energy (per mole of 
H2) consumed by the endothermic electrolysis reaction (TAS) 
and that produced by the irreversible losses within the cell 
(2F ntotal), Which are illustrated in Fig. 6 as a function of current 
density at 1023 K. Note that these results are produced using 
only the electrochemical model in Section 2.1. The cathode 
stream composition, averaged over the cell length, is taken as 
46 mol% H2/54 mol% H20, and employed here together with 
the anode stream composition of 100 mol% O2. As can be 
seen, at 5000 A m7? the irreversible losses are not sufficient 
to provide the entire thermal energy consumed by the reac- 
tion. Consequently, the temperature is decreased as the reaction 
proceeds along the cell. At 7000 Am~*, however, the excess 
heat produced via the irreversible losses causes the temperature 
to rise towards the cell outlet. The curves intersect at around 
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Fig. 6. Thermal energy (per mole of H2) consumed by the reaction (TAS) and 
that produced by the irreversible losses (2F nota) as a function of current density 
at 1023 K. The stream compositions are taken as 46 mol% H2/54 mol% H20 
and 100 mol% Oo. 


5500 A m~*, where the thermal energy consumed by the reac- 
tion is equivalent to that produced by the irreversible losses. Such 
condition corresponds to a thermoneutral operation in which 
the entire enthalpy change in the electrolysis is matched by the 
electrical energy input to the cell. 

In the view of energy efficiency and operating cost, an SOEC 
stack should preferably be driven at current densities below that 
of the thermoneutral operation. Such a condition, which results 
in an endothermic stack operation, allows the opportunity to 
recover waste heat from external processes to partly replace the 
electrical energy consumption of the stack. However, as men- 
tioned previously, the operation of the stack at lower current 
densities requires a larger stack area to support a given H? pro- 
duction rate, involving higher investment cost. The use of heat 
from external sources also involves extra capital cost in the 
form of additional equipment. One point to note at this stage 
is that the curves on Fig. 6 correspond only to the cell opera- 
tion under specified temperature and stream compositions. In an 
actual SOEC, these variables as well as the local current densi- 
ties change along the cell. Therefore, imposing an endothermic 
condition, for example, at the inlet would not necessarily ensure 
the entire cell to operate under the same mode. 

As previously mentioned, the dependence of the cell tem- 
perature distribution on the average current density poses some 
challenges in terms of stack temperature control. Although the 
operating point of an SOEC stack is not expected to change 
dramatically if using grid or nuclear electrical powers as the 
input, such is not the case if the input power varies significantly 
with time. It is thus essential to have an appropriate temperature 
control strategy in place. The implementation of an effective 
control strategy should avoid significant thermal excursions in 
the stack, preventing the fracture of its components and ensuring 
that the operating mode of the stack remains constant. Avoid- 
ing transitions in the operating mode is important as the system 
requirements vary depending on whether the stack is employed 
in an endothermic, thermoneutral or exothermic operation in 
which the average current density is above the thermoneutral 
point. One possible solution in achieving the control require- 
ments would be to modify the SOEC design such that one more 
degree of freedom is available to control the stack temperature. 
As can be seen in [16], in the case of SOFCs, the stack tem- 
perature is often controlled by varying the air flow through the 
cell. Such strategy may potentially be applied in the case of 
an SOEC where the introduction of air flow through the cell 
could both force the O2 product out of the anode channel and 
aid in controlling the temperature. Such air flow would, however, 
cause the stack temperature to decrease and its effect needs to be 
quantified. Different SOEC concepts with the aim of address- 
ing the issue of temperature control are the focus of a future 
publication. 

Figs. 7 and 8 show the contributions of the reversible poten- 
tial and irreversible losses on the cell potential for the conditions 
in Table 3, at the average current densities of 5000 A m~? 
and 7000 A m7, respectively. In both cases, the activation 
overpotentials dominate the irreversible losses while the con- 
centration overpotentials remain insignificant along the cell. For 
the 7000 A m~? case, the thermodynamics predict the reversible 
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Fig. 7. Reversible potential and irreversible losses along the cell for the average 
current density of 5000 A m~? and for the conditions in Table 3. 


potential to decrease as the cell temperature is raised along the 
cell. Fig. 8, however, indicates an increase in the reversible 
potential along the cell due to the change in the cathode stream 
composition. The similarity between the reversible potential dis- 
tributions in Figs. 7 and 8, in spite of the difference in the 
temperature distributions for the two cases, suggests that the 
reversible potential is more strongly coupled with the stream 
composition than with the cell temperature in the range of the 
analysis. Figs. 7 and 8 also show the overall reduction in the irre- 
versible losses along the cell. Although a negative temperature 
gradient is observed in Fig. 5 for the average current density of 
5000 A m~*, Ohmic losses and the activation overpotentials are 
decreased along the cell in Fig. 7. This is due to the rapid reduc- 
tion in the local current density, indicating a significant decrease 
in the Hp production rate along the cell. The local current density 
distributions are shown in Fig. 9 for the average current densities 
of 5000 A m~? and 7000 A m~?. 

The steady state simulations of the full IT-SOEC model for 
the conditions in Table 3 predict cell potentials of 1.28 V and 
1.30V for the average current densities of 5000 A m7? and 
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Fig. 8. Reversible potential and irreversible losses along the cell for the average 
current density of 7000 A m~? and for the conditions in Table 3. 
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Fig. 9. Local current density along the cell for the average current densities of 
5000 A m~? and 7000 A m7”, and for the conditions in Table 3. 


7000 A m~?, respectively. For the 7000 A m~? case in which the 
outlet gas temperature is above the inlet temperature, it may be 
possible to bring the inlet stream up to the required temperature 
entirely through heat recovery from the outlet streams. With 
such a system design, electricity is the only source of energy 
required for the stacks. The predicted 1.30 V equates to an elec- 
tricity consumption of around 250.1 kJ mol~! H3 or3.1 kW h per 
normal m? of H2. This is significantly lower than 4.2 kW h per 
normal m? of H2, which is the electricity consumption of one 
of the most efficient low temperature electrolysis stacks com- 
mercially available today [22]. Neither values include parasitic 
energy consumption within the system. 


4. Conclusions 


A one-dimensional distributed dynamic model of a cathode- 
supported planar IT-SOEC stack has been presented. The model 
has been employed to study both the steady state behaviour 
of an SOEC as well as its electrochemical performance at 
various current densities and temperatures. Electrochemical per- 
formance analysis alone has shown that, for the case where the 
cathode and anode inlet stream compositions are equivalent to 
10 mol% H2/90 mol% H20 and 100 mol% Oz respectively, both 
a decrease in current density and an increase in the cell temper- 
ature result in a reduction of the combined irreversible losses 
and thus a decrease in the electrical potential. Activation over- 
potentials were observed to provide the largest contributions 
to the irreversible losses while the concentration overpotentials 
remained insignificant throughout, in spite of the electrode- 
supported nature of the cell. The steady state performance 
analysis of the cell has shown that for an average current density 
of 7000 A m~? and a cathode inlet stream temperature of 1023 K, 
the cell potential is 1.30 V, which corresponds to an electricity 
consumption of around 3.1 kW h per normal m° of Hp. This is 
significantly smaller than those of low temperature stacks com- 
mercially available today. However, the dependence of the cell 
temperature distribution on the average current density calls for 
a strict temperature control, especially during dynamic opera- 


136 J. Udagawa et al. / Journal of Power Sources 166 (2007) 127-136 


tion. Alternative concepts to accommodate such a feature and 
their dynamic behaviour will be presented in the future. 
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